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ABSTRACT 

We present high spectral resolution (~3 km s _1 ) observations of the V 2 ro-vibrational band of 
H 2 O in the 6.086-6.135 ^irn range toward the massive protostar AFGL 2591 using the Echelon-Cross- 
Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). Ten 
absorption features are detected in total, with seven caused by transitions in the ^2 band of H 2 O, 
two by transitions in the first vibrationally excited 1/2 band of H 2 O, and one by a transition in the 
V 2 band of H^O. Among the detected transitions is the ^2 li,i“0o,o line which probes the lowest 
lying rotational level of para- H 2 O. The stronger transitions appear to be optically thick, but reach 
maximum absorption at a depth of about 25%, suggesting that the background source is only partially 
covered by the absorbing gas, or that the absorption arises within the 6 pm emitting photosphere. 
Assuming a covering fraction of 25%, the H 2 O column density and rotational temperature that best 
fit the observed absorption lines are IV(H20) = (1.3 ± 0.3) x 10 19 cm -2 and T = 640 ± 80 K. 


1. INTRODUCTION 

Water, despite being one of the most abundant species 
in the molecular interstellar medium (ISM), is difficult to 
observe in astrophysical ob jects due to its prevalenc e in 
the Earth’s atmosphere (see lvan Dishocck ct al.ll2013l . for 
a comprehensive review of astronomical water observa¬ 
tions). Ground-based observations of H 2 O have primar¬ 
ily targeted maser emission, most frequently the Jko,k c = 
6i, 6”52,3 transition near 22 GHz that was u tili zed in 
the i nitial detection of interstellar water dCheung et all 
I1969H , although some have also focused on rotational and 
ro-vibrational tra nsitions out of high-lying r otational lev¬ 
els in the mid-IR (iPontomridan et all 1201 (III and near-IR 
dNaiita et all 120001 iC'arr et all 12004 iSalvk et al.1 120081 
llndriolo et alJl2013lf . respectively. These latter observa¬ 
tions are possible because many high-lying rotational lev¬ 
els are not significantly populated in the Earth’s atmo¬ 
sphere, but their scope is limited to astrophysical sources 
with warm (T > 300 K), dense gas. 

Observations of the lowest-lying rotational levels 
of water—those able to probe cold gas—have re¬ 
quired space-based observatories. The Infrared Space 
Observatory -Short Wavelength Spectromete r (ISO- SWS; 
iKessler et al.1119961: Ide Graauw et all I1996D covered the 
V 2 ro-vibrational band (symmetric bending mode) 
of H 2 O centered near 6 pm, and absorption out 
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of the lowest-lying levels of the ortho and para 
nuclear spin modifications (lo.i and 0o,o, respec- 
tivcly) was detected toward sev eral massive protostars 
dBoonman fe van Dishocckll2003f) . Due to the low spec¬ 
tral resolution of the observing configuration though 
(A/AA ~ 1400 using SWS in AOT6 grating mode), 
these lines were significantly blended with absorption 
from other nearby H 2 O lines, making the determina¬ 
tion of level-specific column densities impossible. In¬ 
stead, the entire V 2 band was fit simultaneously assum¬ 
ing a single temperature to determine the total water 
column density, iV(H 2 0). T he Submillimeter Wave As¬ 
tronomy Satellite CSWAS: lMelnick et al.ll2000bh provided 
much higher spectral resolution (< 1 km s~ 1 2 3 4 ) and cov¬ 
ered the li,o~lo.i pure rotational transition of H 2 O at 
557 GHz. This line was observed in both emission and 
absorpti on in multiple sources (e.g. JMclnick et aO2000at 
iSnell et aLll2000f) , demonstrating the ability to probe cold 
water. More recently, the study of low-lying rotational 
levels at high spectral resolution (~ 0.5 km s" 1 ) has 
been facilitated by the Heterodyne Instru ment for the 
Far-Infrared IHIFL lde Graauw et al.ll 2010H on boa rd the 
Herschel Space Observatory (jPilbratt et aLll2010fl . Wa¬ 
ter has been detected in both emission and absorption 
out of levels with E < 200 K in several protostars (e.g., 
Ivan Dishoeck et all 120111 Ivan der Tak et all 1 201 3ll . and 
in absorption o ut of the Op .n and ln.i levels in the molec¬ 
ular ISM (e.g., iSonnentrucker et all 120101 : iFlagev et al.1 
120131 1. Observations that resolve the velocity structure 
of absorption lines are vital to both determining level- 
specific column densities, and understanding the dynam¬ 
ics of the absorbing/emitting regions. This is especially 
important for protostars as such objects contain mul¬ 
tiple dynamical components (e.g., disk, envelope, jets, 
outflows, shocks). 

AFGL 2591 is a region of ongoing high-mass star for¬ 
mation. The radio continuum source VLA 3 (for a de¬ 
scription of sources see iTorrelles et al.l 12014 and refer¬ 
ences therein) is the brightest mid-IR source in the region 
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Figure 1. Bottom: Spectra of AFGL 2591 (black) and Vega (gray) scaled to match in arbitrary intensity units. Water absorption 
features in AFGL 2591 are marked by vertical lines and labeled accordingly. They are blue-shifted by about 40 km s —1 with respect to 
their telluric counterparts due to the Earth’s motion and gas velocity. Top: Spectrum of AFGL 2591 following division by boxcar averaged 
Vega spectrum and normalization. Gaps correspond to regions where the atmosphere is opaque, and poor removal of moderately strong 
atmospheric lines can be seen as spikes to the right of the astrophysical lines (e.g., near 6.089 /mi. 


and drives a bipolar outflow. It was previously observed 
at 5- 7 /im as part of the aforementione d ISO - SWS study, 
from which iBoonman fe van Dishoecd (120031) reported a 
best-fit water column density of -/V(H 2 0) = (3.5 ± 1.5) x 
10 18 cm -2 at T = 450+ 2 gg K, assuming a Doppler line 
width of 5 km s — . Known kinematic components associ¬ 
ated with AFGL 2591 include: (1) the protostellar enve¬ 
lope with systemic velocit y —5.5 km s -1 in the loca l stan¬ 
dard of rest (LSR) frame (Ivan der Tak et al~lll999l f: (2) a 
blue-shifted outflow at —25 km s -1 5. ut.sr 5. —6 km s -1 

(jEmprechtinerer et alJl2012i : lvan der Tak et alll2013H : (3) 

a red-shifted outfl ow at 3 km s -1 < ulsr ^ 15 km s -1 
( Lada et all 1984 ): (4) material entrained in a jet or 
wind, indicated by absorption in the v = 1-0 band 
of 12 CO that extends to —1 96 km s -1 in the form 
of an asymmetric blue wing (Ivan der Tak et ahl I1999D . 
There is also evidence for a rotating disk around 
VLA 3, although gas velocities in the disk match those 
of the envelope (Ivan dcr Tak et all 120061 IWang et al.l 
l2012f) . The interested reader is refe r red to cartoon pic¬ 
tures (Ivan der Tak et all 119991120061 Ivan der Wiel et' al.1 
120111 iSanna et al.l 120121 IWang et alJ~ 2012ft and multi- 
wavelength images (|Johnston et al.ll2013f) of AFGL 2591 
to gain a better understanding of the region. Because 
of the low spectral resolution of the ISO -SWS observa¬ 
tions, it has been impossible to say with certainty which 
component gives rise to the H 2 0 absorption. To do 
so—and to better constrain the water column density 
and rotational temperature—we have targeted multiple 
transitions in the i/ 2 band of H 2 O using the Echelon - 
Cross-Echelle Spectrograph (EXES; iRichter et al.ll20 lDl 
on board the Stra tospheric Observat ory for Infrared As¬ 
tronomy (SOFIA: lYoung et al.ll2012T l. 

SOFIA operates at altitudes above 39,000 ft (11887 m), 


where the precipitable water vapor overburden is rou¬ 
tinely less than 0.02 mm. Under these conditions the 
6 /im region of the Earth’s atmosphere is no longer 
opaque, as it is from the ground. EXES provides high 
spectral resolution (~ 3 km s -1 ) capabilities in the 4.5- 
28.3 /im range, making it well-suited for velocity-resolved 
observations of individual ro-vibrational transitions of 
the i >2 band of H 2 0 in astrophysical sources. We present 
here the first spectrally resolved detections of 10 absorp¬ 
tion lines from transitions in the i/ 2 bands of H 2 0 and 
H 2 8 0, including a detection probing the ground para 
level, 0o,o- This highlights the opportunity to further 
probe water in cold molecular clouds without the need for 
a space-based observatory by utilizing EXES on SOFIA. 

2. OBSERVATIONS AND DATA REDUCTION 

AFGL 2591 VLA 3 was observed using EXES on 
board SOFIA at an altitude of 43000 ft (13106 m) 
on 2014 Apr 10 (UT) as part of instrument commis¬ 
sioning observations. Spectra were acquired in cross- 
dispersed high-resolution mode with a central wavelength 
of 6.1125 /im, using a slit length of 9"9, and a slit width 
of F/9 to provide a resolving power (resolution) of 86,000 
(3.5 km s -1 ), with the resolution element sampled by 8 
pixels. The telescope was nodded after every 27 seconds 
of integration, enabling subtraction of telluric emission 
lines, and the total exposure time for AFGL 2591 was 
1134 s. Prior to observing the target, a calibration se¬ 
quence was taken using the same wavelength setting and 
slit width. This sequence consisted of observations of the 
inter nal blackbody un it set to 260±0.1 K and of blank 
sky (lLacv et al.l I2002T) . These calibration frames were 
used to correct for blaze efficiency, pixel-to-pixel sensi¬ 
tivity variations, and also to provide a first order flux 
calibration using the expected photon counts from the 
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blackbody unit. The bright star Vega was observed dur¬ 
ing a flight leg 2 hr before the science observations at 
the same altitude and air mass as AFGL 2591 for use 
as a telluric standard star. The observing sequence was 
similar to that employed for AFGL 2591. 

Data were proc essed using the Redux pipeline 
(IClarke et al.l l2014f) with the fspextool software 
package—a modifica tion of the Spextool package 
(iCushing et al.l 120041) —which performs source profile 
construction, extraction and background aperture def¬ 
inition, optimal extraction, and wavelength calibration 
for EXES data. The preliminary wavelength scale out¬ 
put from the pipeline was refined by shifting the wave¬ 
length solution so that telluric water absorption lines in 
the AFGL 2591 con tinuum best matched entries in the 
HITRAN database (Rothman et al.l I2013T ). This wave¬ 
length calibration is accurate to within ±0.3 km s -1 . 
Spectra extracted from individual orders of the echcl- 
logram were combined to form a single continuous spec¬ 
trum in the wavelength range 6.085-6.135 pm for both 
AFGL 2591 and Vega as shown in the bottom panel of 
Figure [T] To remove atmospheric features, the spectrum 
of AFGL 2591 was divided by a smoothed (11 point box¬ 
car average) Vega spectrum. The resulting ratioed spec¬ 
trum was then divided by a boxcar average to the con¬ 
tinuum level extrapolated across absorption lines to pro¬ 
duce the normalized spectrum in the top panel of Figure 
H] Figure [2] shows the spectra for each transition shifted 
into the LSR frame following conversion from wavelength 
to line-of-sight velocity. 


3. ANALYSIS 

Absorption lines in Figure [2] were fit using a function 
of the form 


1 = 1 o 


1 ~fc 


1 — exp To exp I — 


(v - ulsr) 5 


2<J?, 


(!) 

where Iq is the continuum level, f c is the fraction of the 
background source covered by absorbing material, and 
the optical depth is assumed to have a Gaussian profile 
with optical depth at line center To, LSR velocity at line 
center ulsr> and velocity dispersion <r„. The V 2 3i,2~ 
22 ,i and U 2 42,2-33 ,\ transitions were fit simultaneously 
(second and fourth panels from the bottom in Figure 
m using a version of this function modified for multiple 
lines. The V 2 53 , 3 - 44,0 and H-^O 1^2 32,i~3i,2 transitions 
(third-from-top and top panels) were fit separately de¬ 
spite their proximity to each other as there is no signifi¬ 
cant blending of the features. Due to interfering features 
near the v 2 li,i-0o,o and V 2 3i,2 _ 3o,3 transitions (bottom 
and third-from-bottom panels; see figure caption for ex¬ 
planation of features) only data at ulsr > — 23 km s^ 1 
were used to constrain fits of these transitions. As the 
absorption features are broad with respect to the instru¬ 
mental spectral resolution, we expect the lines to be re¬ 
solved and convert the optical depth profile fit to a col¬ 
umn density profile via 


dN/dv = r(v (2) 

g u 

under the assumption that the absorption is unsaturated, 
where gi and g u are statistical weights in the lower and 
upper states respectively, A is the spontaneous emission 
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Figure 2. Individual H 2 O absorption lines as functions of LSR 
velocity, arranged in order of increasing lower-state energy (except 
the H^O line). The vertical dashed line shows the systemic ve¬ 
locity of AFGL 2591 deter mined from various mol ecular emission 
lines (ulsr = —5.5 km s~ 1 ; lvan der Tak et al.119991) . Red and blue 
curves show fits with f c = 1 and f c = 0.25, respectively. For opti¬ 
cally thin lines both fits are nearly identical, so only the f c — 1 fits 
are shown. At ulsr, < —25 km s —1 in the v 2 li,i-0o,o spectrum we 
consider the lower continuum level to be caused by poor removal 
of strong atmospheric lines. The arrow at —22 km s —1 marks the 
expected position of the V 2 = 2-1 3i, 3 —2o,2 transition given a line 
center of —11 km s —1 . In the V 2 3i,2-3o, 3 spectrum the arrow at 
—41 km s —1 marks the expected position of the V 2 = 2-1 40,4-31,3 
transition, and we consider the feature from about —32 km s —1 
to —47 km s —1 to be a detection of this line, although it is also 
affected by a sloping baseline due to poor removal of atmospheric 
lines. 


coefficient, and A is the transition wavelength. We then 
integrate over the absorption feature in velocity space to 
determine the column density in the lower state of the 
observed transition. 

Initial fits were made assuming that the absorbing gas 
completely covers the background source (/ c = 1), and 
are shown as red curves in Figure |2j Resulting parame¬ 
ters for each transition are reported in Table [T] Inferred 
column densities are converted to \n{f c N/gi) and plot¬ 
ted versus lower state energy as red squares in Figure [3j 
The black dashed line shows the relationship expected for 
iV(H20) = 3.5 x 10 18 cm -2 in local thermodynamic equi- 
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Table 1 

Observed Water Transitions and Inferred Parameters 


Transition 3, 

Wavelength 

Ei/kt 

91 

A 

^LSR 

(7 v 

TO 

In (fcNi/gi) 


w 

(K) 


(s- 1 ) 

(km s —1 

) (km s _1 ) 









Fit Results for f c 

= 1 


1-0 5 3 ,2-4 4 ,i 

6.0887005 

702.3 

27 

0.35 

- 11.4 

6.0 

0.13T0.03 

35.2±0.3 

1-0 5 2 , 4 - 43,1 

6.0964081 

552.3 

9 

0.73 

- 11.7 

5.5 

0.13T0.03 

35.5T0.3 

2- 1 22,i-2i,2 

6.1009690 

2412.9 

15 

5.92 

- 9.7 

5.5 

0.09T0.03 

32.7T0.4 

H^O 1-0 3 2 ,i-3i, 2 

6.1034870 

248.7 

21 

6.40 

-10.4 

4.7 

0.06T0.03 

31.8T0.5 

1-0 53 , 3 - 44,0 

6.1039868 

702.3 

9 

0.34 

-9.6 

5.5 

0.06T0.03 

35.4±0.6 

1-0 42,2-33,1 

6.1061925 

410.4 

7 

0.53 

-10.5 

7.3 

0.12±0.03 

36.2±0.3 

1-0 3i, 2 - 22,1 

6.1068262 

194.1 

15 

1.12 

-14.2 

8.5 

0.20T0.04 

35.2±0.2 

2- 1 4o,4-3l,3 

6.1131638 

2502.7 

7 

15.8 

-8.4 

4.4 

0.12±0.03 

32.2±0.3 

1-0 3l,2-3o,3 

6.1137707 

196.8 

21 

6.24 

-16.9 

8.5 

0.34±0.04 

34.0T0.1 

1-0 li,i-Oo,o 

6.1163311 

0 

1 

7.46 

-15.7 

8.7 

0.30T0.04 

35.7T0.1 







Fit Results for / c = 

: 0.25 


1-0 53 , 2 - 44,1 

6.0887005 

702.3 

27 

0.35 

-11.4 

5.6 

0.68 jig'll 

35.4t°'i 

1-0 52,4-43,1 

6.0964081 

552.3 

9 

0.73 

-11.8 

5.2 

0-68 jig'll 

35.6±°'i 

2-1 22,i-2i,2 

6.1009690 

2412.9 

15 

5.92 

-9.8 

5.3 

o.44±S;?J 

32.8±°® 

H| 8 0 1-0 3 2 ,i-3i,2 

6.1034870 

248.7 

21 

6.40 

-10.4 

4.5 

0 29~*~ 017 
u - zy -0.15 

3i-8±g;g 

1-0 53 , 3 - 44,0 

6.1039868 

702.3 

9 

0.34 

-9.6 

5.3 

0 25"*” 0 ' 17 
u - zo —0.14 

35.4 jiQg 

1 0 42,2-33,1 

6.1061925 

410.4 

7 

0.53 

-10.3 

6.9 

0 62"*" 0 ’ 25 
u.oz_o 20 

36.4t°'^ 

1-0 3i, 2 - 22,1 

6.1068262 

194.1 

15 

1.12 

-14.2 

7.3 

1 20"*”°' 51 

35.5t°;0 

2- 1 4o,4-3l,3 

6.1131638 

2502.7 

7 

15.8 

-8.4 

4.1 

n 57+ 0 - 24 

u - 0 ' -0.19 

32-4 ±g;| 

1-0 3l,2-3o,3 

6.1137707 

196.8 

21 

6.24 

-16.8 

5.7 

> 5.1 

> 35.0 

1-0 li,i-Oo,o 

6.1163311 

0 

1 

7.46 

-15.6 

6.5 

> 3.0 

> 36.3 


Note. — Columns 2, 3, 4, and 5 give the transition wavelength, A, lower state energy, Ei/k &, lower-state statistical weight, gi, and the spontaneous 
emission coefficient, A, respectively. Columns 6, 7, and 8 give the best-fit parameters from equation Q. with covering fractions of 100% and 25% 
assumed in the top and bottom halves of the table, respectively. Uncertainties in ulsr and a v are estimated to be 0.3 km s —1 at the la level. 
Uncertainties in tq assume that the uncertainty in the absorption depth at line center is equal to the root mean square noise level of 0.03 in the 
continuum of the ratioed spectrum prior to the removal of baseline fluctuations. The Doppler parameter, 6, and line full width at half maximum, 
FWHM, are related to a v via the equations b — a v y /2 and FWHM = a v 2-^/2 ln(2). 

a Transition labels are given as v^—v'^ J' K , K t — J’k"k" "> where a single prime denotes the upper state and double prime denotes the lower state, and 
both the v\ and V 3 vibrational quantum numbers are omitted as they are 0 in all cases. 


librium (LTE) at T = 450 K dBoonman fc van Dishoeckl 
120031 ). Column densities in the 33 , 1 , 43 , 1 , 44 , 1 , and 44,0 
states are in agreement with predictions based on these 
values, but the states with E/ks < 300 K are below pre¬ 
dicted values, while vibrationally excited states are above 
predicted values. No single temperature provides a good 
fit to all nine points, but if we exclude levels with E/ks < 
300 K then we find N(R 2 0) = (3.7 ±0.8) x 10 18 cm -2 in 
LTE at T = 590 ± 50 K (marked by the red dotted line 
in Figure [3]). 

Assuming _ T_ _and_ iV(H20) from 

iBoonman fc van Dishoeckl (|2003t l the Vi 3i,2-3o,3 
transition should be saturated at line center, which 
is not the case. If the background source is partially 
covered by absorbing material, or if the water absorption 
arises within the 6 pm emitting photosphere, then it 
is possible for optically thick, saturated lines to cause 
only a fractional decrease in the continuum level. 
The minimum possible fractional coverage is equal 
to the depth of the strongest absorption line under 
consideration, in our case the v 2 3 i , 2 - 3 o ,3 transition. 
However, the perceived astrophysical absorption from 
this transition is greatly affected by its telluric coun¬ 
terpart, and the line depth is highly dependent on the 
atmospheric division procedure. The same is true for 
the next strongest line, the v 2 li,i-0o,o transition. We 
use both lines in estimating the minimum covering 
fraction, choosing f c = 0.25, but caution that this is 
highly uncertain. Fits using f c = 0.25 are shown as blue 
curves in Figure [21 and resulting parameters are again 


in Table [T] Inferred column densities are plotted as 
open blue diamonds in Figure |3] Scaling by the covering 
fraction demonstrates that column densities inferred 
from optically thick transitions increase much more 
than those inferred from optically thin transitions as f c 
decreases, but at f c = 0.25 the saturated v 2 31,2-30,3 
and u 2 li.i-0o,o transitions only provide lower limits on 
column densities. A fit to the 7 unsaturated transitions 
is shown by the blue dash-dot line in Figure [H and 
corresponds to iV(H20) = (1.3 ± 0.3) x 10 19 cm -2 in 
LTE at T = 640 ± 80 K. 

4. DISCUSSION 

The pure rotational li,i-0o,o transition of H 2 O at 
1113.343 GHz (269.3 pm) is se en in both emission and 
absorption toward AFGL 2591 (Iv an der Tak et al Jl2013t 

iKazmierczak-Barthel et (±112014 iChoi et al.ll2015l) . The 

emission component arises in the protostcllar envelope, 
and the absorption components in the blue-shifted out¬ 
flow and foreground gas. This absorption probes the 
same quantum state as the u 2 li,i-0o,o transition, but 
a direct comparison is hindered by several effects. The 
emitting regions at 6 pm and 1.1 THz (~270 pm) are 
likely different, so that gas probed by one transition may 
not necessarily be probed by the other. Our EXES ob¬ 
servations used a 9"9 by 1"9 slit centered on VLA 3, 
while HIFI observations at 1113 GHz have a roughly cir¬ 
cular beam with FWHM ~ 19", again meaning that 
different regions are being probed. Finally, emission 
from the {v\v 2 v^) Jic o ,if c =(000) lyi state (53.4 K above 
ground) is extremely strong and interferes with the ab- 
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Figure 3. Rotation diagram for the 9 detected H 2 O t ransitions. The black dashed line is for N( H 2 O) = 3.5 x 10 18 cm -2 and LTE at 
450 K as reported by IBoonman &; van Dishoeckl (120031 ). Red squares and open blue diamonds denote covering fractions of 1 and 0.25, 
respectively. Lower limits are indicated by upward-pointing arrows. Note that the y-axis shows ln(/ c 7V( JK a K c )/9K a k c ), (i.e., is scaled 
to the average column density in front of the entire background source). For optically thin absorption lines the value does not change 
significantly with covering fraction, but for optically thick lines (e.g., V 2 li,i - 0o,o) the value increases substantially as f c approaches the 
absorption depth. The red dotted line marks T = 590 K and N(H.20) = 3.7 x 10 18 cm -2 for f c = 1, and is the best fit to levels with 
E/ks > 300 K. The blue dash-dot line marks T = 640 K and iV(H 2 0) = 1.3 x 10 19 cm -2 for f c = 0.25, and is the best fit to the 7 
unsaturated transitions. 


sorption features at 1113 GHz, whereas emission from 
the (010) I 14 state (2352 K above ground) is not ob¬ 
served, leaving the 6.1163311 /.tm line unobscured. All 
of these effects must be considered when comparing any 
results from HIFI and EXES. 

As described in Section [Q different components of 
AFGL 2591 are distinguished by line-of-sight velocities, 
line profiles, and observed molecules. The H 2 O absorp¬ 
tion lines are centered at roughly —16 km s _1 with 
FWHM ~ 20 km s -1 for the three lowest-lying levels 
and shift to —11 km s -1 with FWHM ~ 13 km s -1 for 
the higher-lying levels, best matching the blue-shifted 
outflow. These line profiles are most similar to ab¬ 
sorption seen in HCN and C 2 H 2 dKnezl [20061 observed 
at 13 (im) as well as specific absorption components 
of 13 CO and vibrational ly exc ited 12 CO dMitchell et al.1 
119891 Ivan der Tak et all 119991 observed at 4.7 /zm), all 
of which are thought to arise in hot, dense gas close 
to the central protostar. Several 22 GHz H 2 O masers 
are observed throughout the region, with many concen¬ 
trated in the w alls of the blue-shifted outflow associ¬ 
ated with VLA 3 dTrinidad et al.ll2003l : lSanna et al.ll 2012 l ; 
iTorrelles et al.ll2014f) . Although maser velocities are pri¬ 

marily at ^lsr < —18 km s , absorption by gas in small, 
shocked clumps could explain the small covering fraction 


required by our analysis. Given the above, we posit that 
the H 2 0 absorption observed with EXES arises in hot, 
dense gas at the base of the blue-shifted outflow. 

Water emission attributed to the outflow com¬ 
ponent is observed with H IFI in multiple t ransi¬ 
tions as a broad feature (Ivan der Tak et al.1 120131: 
iKazmierczak-Barthel et ah! 12014 iChoi et a.l.l 1 20 loll , but 

is centered closer to the systemic velocity of the en¬ 
velope than the absorptio n we see. Analy sis of the 
H 2 0 outflow emission by iGhoi et ahl (20151 1 indicates 
a temperature of T ~ 70-90 K and column density 
of y(H-;>0) ~ 4 x 10 13 cm -2 , similar to results found 
by iKarska et al.l (I2014H who analyzed unresolved H 2 0 
absorption observed with PACS on Herschel , finding 
T = 160 ± 130 K and iV(H 2 0) ~ 10 14 cm" 2 . These 
values are significantly below what we find, suggesting 
that EXES and HIFI/PACS observations trace different 
components. Between the broad and narrow water emis¬ 
sion observed with HIFI, the broad absorption observed 
with EXES, and the abundant maser spots, it is evident 
that the water-containing gas in AFGL 2591 is both spa¬ 
tially and kinematically complex. Interpreting the var¬ 
ious observations of H 2 0 in AFGL 2591 in unison will 
require utilizing a physical model that includes radiative 
transfer accounting for absorption, stimulated and spon- 
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taneous emission, and collisional (de-)excitation, as well 
as kinematic and geometric effects. 

5. SUMMARY 

We have detected ten absorption features arising from 
warm gas in AFGL 2591 caused by ro-vibrational tran¬ 
sitions of water, including seven from the V 2 band of 
H 2 O, two from the vibrationally excited band of H 2 O, 
and one from the band of H^O. Among the de¬ 
tected absorption lines is the li.i~0o,o transition at 
6.1163311 /im, which probes the ground state of para- 
H 2 O. Relative strengths of absorption features are sug¬ 
gestive of a covering fraction less than 1 (or absorption 
arising within the 6 emitting photosphere), with a 
limit of f c > 0.25 set by the depth of the strongest ab¬ 
sorption features. Analysis of the level populations as¬ 
suming f c = 0.25 results in V(H 2 0) = (1.3 ± 0.3) x 
10 19 cm -2 for LTE at T = 640 ± 80 K. Line profiles best 
match the blue-shifted outflow component, and we as¬ 
cribe the absorption to hot, dense gas at the base of the 
outflow. The temperature and column density inferred 
by our ana l ysis a re much larger than those reported by 
iChoi et all (|2015H for the outflow component observed 
in H 2 O emission by HIFI. Uncertainty in whether the 
EXES and HIFI observations are probing the same gas 
makes the combined interpretation of both datasets dif¬ 
ficult, and a physical model of AFGL 2591 that includes 
radiative transfer will be necessary for such an analy¬ 
sis. Clearly though, observations of the v 2 ro-vibrational 
band of H 2 O at high spectral resolution—observations 
uniquely achievable with EXES on SOFIA—add impor¬ 
tant information for interpreting this region and others. 

MJR and CND acknowledge Collaborative Agreement 
NNX13AI85A between UCD and NASA Ames for its 
support and support of EXES development. Many 
thanks to the anonymous referee. 
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